A dynamical approach to the spatiotemporal aspects of the Portevin-Le ChateUer effect: Chaos, 
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The analysis of experimental time series obtained from single and poly-crystals subjected to a constant strain 
rate tests report an intriguing dynamical crossover from a low dimensional chaotic state at medium strain rates 
to an infinite dimensional power law state of stress drops at high strain rates. We present results of an extensive 
study of all aspects of the PLC effect within the context a recent model that reproduces this crossover. We 
characterize the dynamics of this crossover by studying the distribution of the Lyapunov exponents as a function 
of the strain rate, with special attention to system size effects. The distribution of the exponents changes from 
a small set of positive exponents in the chaotic regime to a dense set of null exponents in the scaling regime. 
As the latter is similar to the result in the GOY shell model for turbulence, we compare the results of our model 
with that of the GOY model. Interestingly, the null exponents in our model themselves obey a power law. The 
study is complimented by visualizing the configuration of dislocations through the slow manifold analysis. This 
shows that while a large proportion of dislocations are in the pinned state in the chaotic regime, most of them 
are pushed to the threshold of unpinning in the scaling regime, thus providing an insight into the mechanism of 
crossover. We also show that this model qualitatively reproduces the different types of deformation bands seen 
in experiments. At high strain rates where propagating bands are seen, the model equations can be reduced to 
the Fisher-Kolmogorov equation for propagative fronts. Marginal stability analysis shows that the velocity of 
the propagating of the bands varies linearly with the strain rate and inversely with the dislocation density. These 
results are consistent with the known experimental results. We also discuss the connection between the nature 
of band types and the dynamics in the respective regimes. The analysis demonstrates that this simple dynamical 
model captures the complex spatio-temporal features of the PLC effect. 

PACS numbers: 05.65.+b, 05.45.Ac, 62.20.Fe, 05.9().+m 



I. INTRODUCTION 



Plastic deformation is a highly dissipative irreversible 
nonequilibrium process where nonlinearities play a funda- 
mental role. Under normal conditions, one finds homoge- 
neous deformation. However, under suitable conditions of de- 
formation, different types of spatial and temporal patterns are 
observed. At a microscopic level, these patterns arise due to 
the collective behavior of dislocations. These can be broadly 
classified by the associated time and length scales. For exam- 
ple, the persistent slip bands observed in cyclic deformation is 
an exarnple of nearly permanent pattern lasting over long time 
scales On the other hand, a type of propagative bands 
referred to as the Liiders bands observed in uniaxial tension 
tests are characterized by short time scale Ql . Yet another and 
even more complex spatio-temporal patterns are observed dur- 
ing tension tests of dilute metallic alloys in a certain range of 
strain rates and temperatures. This phenomenon has come to 
be known as the Portevin-Le Chatelier (PLC) effect |2]. Here 
a uniform deformation mode becomes unstable leading to a 
spatially and temporally inhomogeneous state. The instability 
manifests itself in the form of serrations on the stress-strain 
curves of the sample (3|3|. Each stress drop is generally asso- 
ciated with the nucleation and often the propagation of a band 
of localized plastic deformation. In poly-crystals, these bands 
and the associated serrations are classified into three generic 
types. On increasing strain rate or decreasing the temperature, 
one first finds the type C band, identified with randomly nucle- 
ated static bands with large characteristic stress drops on the 
stress- strain curve. The serrations are quite regular Then the 



type B 'hopping' bands are seen. The serrations are more ir- 
regular with amplitudes that are smaller than that for the type 
C. The bands formed are still localized and static in nature, 
but forming ahead of the previous band in a spatially corre- 
lated way giving the visual impression of a hopping propaga- 
tion. Finally, one observes the continuously propagating type 
A bands associated with small stress drops. (In single crys- 
tals such a clear classification does not exist.) These different 
types of PLC bands are believed to represent distinct corre- 
lated states of dislocations in the bands. 

The well accepted classical explanation of the PLC effect 
is via the dynamic strain aging (DSA) concept first introduced 
by Cottrell |4] and later extended by others |1, 5, 6, 7|. In 
the Cotrell's picture, the dynamic strain aging refers to the 
interaction of mobile dislocations with the diffusing solute 
atoms. At low strain rates (or high temperatures) the average 
velocity of dislocations is low and there is sufficient time 
for the solute atoms to diffuse to the dislocations and pin 
them (called as aging). Thus, longer the dislocations are 
arrested, larger will be the stress required to unpin them. 
When these dislocations are unpinned, they move at large 
speeds till they are arrested again. At high strain rates (or 
low temperatures), the time available for solute atoms to 
diffuse to the dislocations decreases and hence the stress 
required to unpin them decreases. Thus, in a range of strain 
rates and temperatures where these two time scales are of 
the same order of magnitude, the PLC instability manifests. 
The competition between the slow rate of pinning and 
sudden unpinning of the dislocations, at the macroscopic 
level translates into a negative strain rate sensitivity (SRS) 
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of the flow stress as a function of strain rate which is the 
basic instability mechanism used in most phenomenological 
models Ql|2l- Slow-fast dynamics and the negative flow rate 
characteristic is common to many stick-slip systems such as 
frictional sliding [8.1, fault dynamics [9] and peeling of an 
adhesive tape llflll and charge density waves 

There are two different types of challenges in dealing with 
the PLC effect. First, understanding the collective behavior of 
dislocations which has been slow largely due to lack of tech- 
niques for describing the cooperative behavior of dislocations. 
Second, the PLC effect involves collective modes of dislo- 
cations where both fast and slow times scales play equally 
important role which requires specific techniques of nonlin- 
ear dynamics (as we shall see). Further, these time scales 
themselves evolve as a function of strain rate and tempera- 
ture which in turn leads to different types of serrations. At 
low strain rate, the existence of both fast ( time scales over 
which stress drops occur) and slow time scales ( loading time 
scales) are clearly displayed in the stress -strain curves. How- 
ever, at high strain rate, as internal (plastic) relaxation is not 
complete, and a clear demarkation of time scales becomes dif- 
ficult. This along with the corresponding length scales ( band 
widths), which also evolve, points to an extremely complex 
underlying dynamics. 

The inherent nonlinearity and the presence of multiple time 
scales necessitates the use of tools and concepts of nonlin- 
ear dynamics for a proper understanding of this phenomenon. 
Early theories based on DSA do not deal with the temporal as- 
pect 1 5, 6, 7] and thus are unsuitable for analyzing the dynam- 
ical aspects of the PLC effect. The first dynamical approach 
was undertaken in early 80s by Ananthakrishna and cowork- 
ers [12.1, which by its very nature affords a natural basis for 
the description of the time dependent aspects of the PLC ef- 
fect. Further, it also allows for explicit inclusion and interplay 
of different time scales inherent in the dynamics of disloca- 
tions. The original model which attempts to address the time 
dependence of the phenomenon uses three types of disloca- 
tion densities assumed to represent the collective degrees of 
freedom of dislocations |12]. Despite the simplicity of the 
model, many generic features of the PLC effect such as the 
existence of a window of strain rates and temperatures within 
which it occurs, etc., were correctly reproduced. More impor- 
tantly, the negative SRS was shown to emerge naturally in the 
model, as a result of nonlinear interaction of the participating 
defects fllfll. 

Due to the dynamical nature of the model, one prediction 
that is unique to this model is the existence of the chaotic 
stress drops in a certain range of temperatures and strain rates 
fl4ll . This triggered a series of experiments to verify this pre- 
diction. The method followed was to analyze the stress-time 
series using dynamical methods Ill7lll8ll . Apart from 

confirming the chaotic nature of stress drops in a window of 
strain rates, these attempts have shown that a wealth of dy- 
namical information can be extracted from the stress-time se- 
ries obtained during the PLC effect 1 15, 16]. Indeed, the num- 
ber of degrees of freedom estimated from the experimental 
time series turn out to be same as in the model offering justi- 



fication for ignoring spatial degrees of freedom. Subsequent 
efforts to extend this analysis to the time series obtained over 
a range of strain rates showed an intriguing crossover from a 
chaotic state at low and medium strain rates to a power law 
state at high strain rates [19, 20 ] . As the crossover is observed 
in both single and polycrystals, it appears to be insensitive to 
the microstructure. However, chaotic state is dynamically a 
distinct state from the power law state as the former involves a 
small number of degrees of freedom characterized by the self- 
similarity of the attractor and sensitivity to initial conditions 
1 18] while latter is an infinite dimensional state reminiscent 
of self-organized criticality (SOC) ]21, Due to this 

basic difference in the nature of the dynamics, most systems 
exhibit either of these states. More importantly, these studies 
also demonstrate that the nature of the dynamics in a given 
strain regime is correlated with the nature of band type. The 
chaotic state has been identified with the type B bands and the 
scaling regime at high strain rate with the propagating type 
A bands ]20]. These authors also make a connection between 
the transition in the nature of serration between the type B and 
type A bands regime of strain rates with the Anderson's transi- 
tion in condensed matter physics. Indeed, recently the spatio- 
temporal features of the PLC effect have attracted attention 
from physicists also L24J . Thus, it appears that the PLC ef- 
fect is a storehouse of many paradigms in condensed matter 
physics. Understanding these connections between dynamics 
and general features of the PLC effect would give insight into 
the rich physics. As the above studies underscore the impor- 
tance of nonlinearity, it demands a dynamical approach to the 
PLC effect. 

The dynamics of the crossover as a function of strain rate 
is unusual in a number of ways. First, the PLC effect is one 
of the two rare instances where such an intriguing crossover 
phenomenon is seen, the other being in the hydrodynamic tur- 
bulence f2^. Second, the power law, both in the PLC ef- 
fect and turbulence, arises at high drive rates ]26, 27]. Thus, 
it would be interesting to examine the similarity and differ- 
ences with hydrodynamic turbulence by comparing results of 
the Lyapunov spectrum of the model for the PLC effect and 
GOY shell model of turbulence {tA 12811 . Another motivation 
is that such a study helps us to compare the nature of the Lya- 
punov spectrum with the conventional SOC systems seen at 
low drives ( such as those in earthquakes ]29], acoustic emis- 
sion during volcanic activity ]30], Barkhausen noise ]31j). ( 
For lack of anything better, we shall reserve SOC for power 
law situations at low drives.) Finally, as different types of 
bands are a characteristic feature of the PLC effect, we inves- 
tigate the connection between spatial aspects and the nature of 
the dynamics. 

The fully dynamical nature of the model and it prediction of 
chaotic stress drops at intermediate strain rates as found in ex- 
periments, makes it most suitable for studying this crossover 
by including spatial degrees of freedom. The paper reports a 
detailed investigation of all these issues ( some of which has 
been reported in brief earlier fs^jIsS*]) within the context of an 
extension of Ananthakrishna's model for the PLC effect. Par- 
ticular attention will be paid to study the system size effects 
during the crossover. 
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Section II, briefly introduces the dynamical model and its 
extension to include spatial degrees of freedom. Section III 
contains the numerical procedure used. In Section IV, we in- 
troduces the background material used for the study. Section 
V contains a comparison of the results of analysis of experi- 
mental time series with that of the model. Section VI, contains 
all the major results on the dynamics of crossover including 
the evolution of the Lyapunov spectrum with special attention 
to study the system size effects as a function of the strain rate 
along with the analysis of the distribution of null exponents 
in the power law regime of stress drops. The section also in- 
cludes a comparison of the results of the model with that of 
the GOY model for turbulence followed by the slow manifold 
method of visual realization of dislocation configurations. Fi- 
nally, in section VII we discuss both analytical and numerical 
results on the nature of dislocation bands. We conclude the 
paper with a few general comments. 



II. THE ANANTHAKRSHANA'S MODEL 

In the dynamical model due to Ananthakrishna and cowork- 
ers fT2il . the well separated time scales mentioned in the DSA 
are mimicked by three types of dislocations, namely, the fast 
mobile, immobile and the 'decorated' Cottrell type disloca- 
tions. The basic idea of the model is that all the qualitative 
features of the PLC effect emerge from the nonlinear inter- 
action of these few dislocation populations, assumed to rep- 
resent the collective degrees of freedom of the system. As 
the model has been studied in detail by our group and others 
including an extension to the case of fatigue |34, 35, 3^, fol- 
lowing the notation in Ref. iTsll . we shall briefly outline the 
model in the scaled variables. In our model, a natural basis for 
including the spatial coupling is through the cross-slip mech- 
anism proposed earlier 1 1 ] with an important difference (see 
below). The model consists of densities of mobile, immobile, 
and Cottrell's type dislocations denoted by p,„{x,t), p,„,(x,f) 
and Pc{x,t) respectively, in the scaled form. The evolution 
equations are: 



dpm 
dt 



dt 

dpc 

dt 



-bop,,, - p„,p;„, + p,„, - ap,„ + (i?'"ffP„ 

D_dH^^^^ix)p^ 
Pi,„ dx^ 



iPm-pmPm-pim+aPc), 



c(Pm - Pc) 



(1) 
(2) 
(3) 



The model includes the following dislocation mechanisms: 
immobilization of two mobile dislocations due to the forma- 
tion of locks (bopl,), the annihilation of a mobile dislocation 
with an immobile one {pmPim), the remobilisation of the im- 
mobile dislocation due to stress or thermal activation (p;,„). It 
also includes the immobilisation of mobile dislocations due to 
solute atoms (op,,,). Once a mobile dislocation starts acquiring 
solute atoms we regard it as the Cottrell's type dislocation p^. 



As they progressively acquire more solute atoms, they even- 
tually stop, then they are considered as immobile dislocations 
pi,n- Alternately, the aggregation of solute atoms can be re- 
garded as the definition of p^, ie., p^ — f'_„dt'p,„{t')K{t — t'), 
where K{t) is an appropriate kernel. For the sake of simplicity, 
this kernel is modelled through a single time scale, K{t) = e^'^' 
. The convoluted nature of the integral physically implies that 
the mobile dislocations to which solute atoms aggregate ear- 
lier will be aged more than those which acquire solute atoms 
later (see ref. 1 13]). The fifth term in Eqn.(l) represents the 
rate of multiplication of dislocations due to cross-slip. This 
depends on the velocity of the mobile dislocations taken to be 

Vm{<^) = ¥!^ffy where (^eff = {^^hp]l'^) is the scaled effec- 
tive stress, (|) the scaled stress, m the velocity exponent and h 
a work hardening parameter. 

The nature of the spatial coupling in the PLC effect has 
been a matter of much debate 1 1]. Several mechanisms have 
been suggested as a source of spatial coupling, such as com- 
patibiUty stresses between the slipped and the un-slipped re- 
gions, long range interactions, and triaxiality of stresses (l|. 
Within the scope of our model, cross-slip is a natural source of 
spatial coupling, as dislocations generated due to cross slip at 
a point spread over to the neighboring elements. Let Ax be an 
elementary length. Then, the flux <f>(x) flowing from jciAx 
and out of x is given by 



4>(x) + |[4>(x- 



-Ax)-2^{x)+^{x-Ax)]). 



(4) 



where (^{x) — Pm(jc)Vm(-x) and p is the probability of cross- 
slip spreading into neighboring elements. Expanding <I>(jf± 
Ax) and keeping the leading terms, we get 



(5) 



We further note that cross-slip spreads only into regions of 
minimum back stress. Here, we consider the back stress is 
taken to result from the immobile dislocation density ahead 
of it. Thus, we use Ax^ =< Ax^ >= r^pJi,^, where < . . . > 
refers to the ensemble average and is an elementary (di- 
mensionless) length. Finally, a, bo and c are the scaled rate 
constants referring, respectively, to the concentration of solute 
atoms slowing down the mobile dislocations, the thermal and 
athermal reactivation of immobile dislocations, and the rate at 
which the solute atoms are gathering around the mobile dis- 
locations. We note here that the order of magnitudes of the 
constants have been identified in Ref fl^ Issl l36ll . These 
equations are coupled to the machine equation 



dt ^ I 



Pm{x,t)<^effix,t)dx], 



(6) 



where e is the scaled applied strain rate, d the scaled effective 
modulus of the machine and the sample, and / the dimension- 
less length of the sample. (We reserve for the unsealed 
strain rate.) We also note here that there is a feed back mecha- 
nism between Eq. |6]and Eq.(l). The machine equation which 
determines the stress depends on the the difference between 
the applied strain rate and average plastic strain rate generated 
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in the sample. Thus, the nature of internal relaxation can influ- 
ence stress generated in the sample which in turn determines 
the dislocation multiplication in Eq. (1). This type of global 
coupling (Eq. |6|l is common to many other situations for in- 
stance in the nonlinear transport properties of charge density 
waves ( in blue bronze for example) 1 1 1]. We shall make some 
comments on this later. 



III. NUMERICAL SOLUTION OF THE MODEL 

We first note that the spatial dependence of p,„, and arises 
only through that of p,„. We solve the above set of equations 
by discretizing the specimen length into equal parts. Then, 
Pm{j,t), Pim{j,t), Pc{j,t), j = ^L^N, and (|)(f) are solved. 
The widely differing time scales Il3l l32. 37] calls for appro- 
priate care in the numerical solutions. We use a variable step 
fourth order Runge-Kutta scheme with an accuracy of 10^^ 
for all the four variables. The spatial derivative in p„j is ap- 
proximated by its central difference. The initial values of the 
dislocation densities are so chosen that they mimic the val- 
ues in real samples. They are uniformly distributed with a 
Gaussian spread along the sample. However, for most calcu- 
lations, we have used the steady state values for the variables 
as the long term evolution does not depend on the initial val- 
ues. As for the boundary conditions, we note that the sample 
is strained at the grips. This means that there is a high density 
of immobile dislocations at the ends of the sample. We simu- 
late this by employing two orders of magnitude higher values 
for Pi,„{j,t) at the end points j = 1, and than the rest of the 
sample. Further, as bands cannot propagate into the grips, we 
use Pm(;,f) = Pc{j,t) at 7 1 and A^. 

For the original model (D = 0), it has been shown that the 
fixed point of the system of equations becomes unstable in a 
certain range of parameter values. In particular, as a function 
of the applied strain rate, the PLC state is reached through 
a Hopf bifurcation and is terminated by a reverse Hopf bi- 
furcation ( with the other parameters kept in the instability 
domain). This feature is retained with the addition of the spa- 
tial degrees of freedom except that the number of complex 
conjugate roots are 2N, the negative ones are N and one zero 
exponent. We find that the instability domain in e increases 
when the values of the other parameters a,bo,c,d, m are taken 
as in Ref. |37]. This is due to fact that the range of e de- 
pends on the value of D due to the global coupling in Eq. 
|6l ( The domain converges quickly as a function of A^.) The 
boundary of e is approximately in the range 10 to 1000 for 
a = 0.8,^0 = 0.0005, c = 0.08, c/ = 0.00006, m = 3.0, /i = 
with D — 0.5, beyond which a uniform steady state exits. A 
set of four eigen values are shown in Fig. ^ The numerical 
results reported in the present work are for the above values. 
However, the results hold true for a wide range of values of 
other parameters in the instability domain including a range 
of values of D. Various system sizes are used depending upon 
the property studied, but are generally in the range A^=100 to 
3333. A sequence of values of N are used wherever conver- 
gence of the properties are investigated. 




FIG. 1: Eigen value spectrum of the fixed point for the model, co,- 
and CO; refer to the real and imaginary parts of the eigen value. 

IV. METHODOLOGY 

As our approach is fully dynamical and keeping in view the 
materials science community, we collect here a few defini- 
tions and provide some details of the methodology used in the 
analysis. Characterizing the dynamics of the model equations 
is carried out by studying the Lyapunov spectrum. The num- 
ber of Lyapunov exponents M for a given N is M — 3N + 1 . 
We shall also use two other well know invariants namely the 

yj y^. 

Kaplan- Yorke dimension Dky — j + \x^l^ \ ^ where / is such 

that L/=i^i > 7^/1=1 '^i < and the Kolmogorov entropy 
H = T.f^i A-i such that Xp > and A-^+i < 0. One important 
issue that is relevant to systems with many degrees of free- 
dom is the existence of a limiting density for the Lyapunov 
spectrum as the system size is increased. This requires that 
we should ascertain if Xj versus x — j/U' converges to a well- 
defined asymptotic density function A(x) with x e [0,1]. ( See 
Ref- L27J.) We address this issue by calculating the spectrum 
for various system sizes A^ = 100 to 3333 which covers ap- 
proximately two orders in M. In particular, such a study will 
be useful in comparing the results of our model with the GOY 
shell model for turbulence ll28il in the power law regime of 
stress drops. Then, one expects that j /Dky converges to well 
defined density function. Following Ref. |28], we use j /Dky 
verses an appropriately scaled quantity XjDky /H. This quan- 
tity is expected to converge to f(kjDKY/H). ( We note here 
that the distribution function is proportional to the negative 
derivative of /.) The nature of the converged Lyapunov den- 
sity function /{XjDky /H) as a function of the drive parameter 
e can be used to quantify the changes the dynamics during the 
crossover. 

As stated earlier, a proper description of the PLC effect 
requires a description of both the slow and fast time scales 
which in turn requires special techniques in nonlinear dynam- 
ics. These two time scales are transparent in the model equa- 
tions where Eq. ([0 represents a fast dynamics compared to 
the rest ( both Eq. (|2j and (|6} are slow while (|3} falls in be- 
tween). Such a system can be studied by eliminating the fast 
variable thereby allowing a reduction in the dimensionality of 
the system |38|. To illustrate this consider 

A« = f{x,y,^') (7) 
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y = 8{x,y,ij) (8) 



where /j is small parameter and x e M'' and y G K''. The main 
feature of such systems is that x evolves much faster than y 
unless f{x,y,/i) is small. In the vicinity of the slow manifold 
defined by f{x,y,/j) — 0, the dynamics is characterized by the 
evolution of the slow variable y. Thus, there is a reduction in 
the dimensionality of the system. On the other hand, if one 
is interested in the fast sub-system, using a scaled time T = 
f /jU, we get the corresponding fast variable x defined by Eq. 
where the slow variables y act as parameters (obtained from 
Eq. |8}. This subspace is clearly the complimentary subspace 
of the slow manifold. We shall use these two subspaces for the 
visualization of dislocation configurations in the high strain 
rate power law regime and obtain the band velocity at high 
strain rates respectively. 




5000 6000 2100 2600 
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FIG. 2: (a) & (b) Experimental stress-time series: (a) chaotic state 
at strain rates = 1.7 x lO^^i'^' and (b) power law state at ta = 
8.3 X lO^^i'^'. (c) & (d) Stress-time series from the model at (c) 
e= 120 (d)e = 280. 



The analysis of the experimental stress-time series is car- 
ried out by estimating both correlation dimension v and the 
Lyapunov spectrum. These methods involve embedding the 
scalar time series in a higher dimensional space using time- 
delay technique fisllzsll . Givenatime series {Oj\j = 1, ..,M}, 

one first constructs vectors ^/ = (o,- , a,_x , 2x , • ■ • , Om- (i/- 1 )x ) 
in a c/ dimensional space. The assumption here is that the 
actual dynamics can be unfolded by embedding the time se- 
ries in a higher dimensional space in which the original at- 
tractor resides. (In addition, surrogate data analysis was 
also carried out in |19].) Then, a quantitative estimate of 
the self similarity of the attractor, namely the correlation di- 
mension, V, can be obtained by calculating the integral l^9il 
C{r) ^ ^L^ir- iC' -1/1) ^ r", where Nt is the total num- 
ber of points in the sum. Correlation dimension also pro- 
vides a lower bound for the number degrees of freedom re- 
quired for a dynamical description of the system which is 
given by the minimum integer larger than v + 1 |40|. The 
geometrical interpretation of these degrees of freedom is that 
they correspond to the subspace to which the trajectories are 
confined. The dimension of this subspace can be obtained 
directly by using singular value decomposition (SVD) |41]. 
This method is often used for filtering noise component su- 
perposed on the time series. However, another use of SVD in 
the present context is that it is useful for the visualization of 
the strange attractor. (This method has been applied to the 
PLC time series earlier |16].) The method involves setting 
up the m X d trajectory matrix T defined by (^1,^2, •••,^m) 
where m — M — {d ~ I) x x. The eigen values of the ma- 
trix are obtained using the standard method of decomposition 
T = UWV'^, where U is m X d orthogonal matrix, \ is adx d 
unitary matrix and W is the matrix of eigen values of the co- 
variance matrix of T which are all nonnegative. The eigen val- 
ues usually decrease rapidly saturating to a level below which 
the changes are minimal. Then the dimension of the attractor 
is taken to be that corresponding to number at which the eigen 
values saturate. 
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FIG. 3: Singular value spectrum of the experimental time series 
shown in Fig.|2t 



V. COMPARISON WITH EXPERIMENTS 

To make the motivation clear, we begin by briefly recall- 
ing the relevant experimental results on the crossover phe- 
nomenon and then compare them with those from the model. 
The simplest feature to compare is the nature of serrations in 
the respective regimes of strain rate. Plots of two experimental 
stress-strain curves from CuAl single crystals corresponding 
to the chaotic and power law regimes of applied strain rates 
are shown in Fig. |2]a,b. The stress-time series in the interme- 
diate and high strain rate regimes from the model are shown in 
Fig-Efc'd. The similarity of the experimental time series with 
that of the model in the respective regimes are clear 

The analysis of the stress-time series given in Fig. |2]has 
been reported in Ref. The correlation dimension was 

found to be V = 2.3. Then, the number of degrees of free- 
dom required for the description of the dynamics of the system 
given by the minimum integer larger than v + 1 |40] is seen to 
be four, consistent with that used in the original model. As an 
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FIG. 4: (a) Reconstructed experimental attractor from the time series 
shown in Fig.|2t. (b)Attractor from the model for A' = 100, j = 50. 



independent check to obtain the number degrees of freedom, 
as also for the visualization of the experimental attractor, we 
have caiTied out the singular value decomposition of this time 
series. The normalized eigen values ( with respect the largest) 
is shown in Fig. |3] It is clear that the relative strength of the 
fourth eigen value drops more than two orders of magnitude 
compared to the first and changes very little beyond the fourth. 
Thus, we estimate the dimension of the experimental attractor 
to be four which is also consistent with that obtained from the 
correlation dimension. ( For time series from model systems, 
one usually finds a floor level below which the eigen values 
saturate. This is taken as the dimension of the actual attractor. 
However, in real situations, as in the present case, the eigen 
values do not saturate due to the presence of noise.) Then, for 
the visualization of the experimental attractor, we can use the 
dominant eigen values to reconstruct the nature of the attrac- 
tor. Using the first three principal directions of the subspace 
C,;/ = 1 to 3, we have reconstructed the experimental attrac- 
tor in the space of specifically chosen directions Ci ~C2,C3 
and Ci to permit comparison with the attractor obtained from 
the model. This is shown in Fig. ^ a for the experimental 
time series ate„ = 1.7 x lO^^i^'. This can be compared with 
the strange attractor obtained from the model in the space of 
Pm,Pim and P( (at an arbitrary spatial location, here j ~ 50 
and = 100) shown in Fig.l^b for 8 = 120 corresponding to 
the mid chaotic region (see below). Note the similarity with 
the experimental attractor particularly about the linear portion 
in the phase space (Fig.|4]a). This direction can be identified 
with the loading direction in Fig.|2]a. Note that the identifica- 
tion of the loading direction is consistent with the absence of 
growth of p,„. 




FIG. 5: Distributions of the stress drops from the model (o), from 
experiments (•) for A' = 1000 and e = 280. Solid line is a guide to 
the eye. 



In contrast to the experimental time series at low and 
medium strain rates, for the time series at the highest strain 
rate ( Fig. Elb), we neither find a positive Lyapunov exponent 
and nor a conv erge d value of the correlation dimension as was 
shown in Ref. lll9ll . However, the distribution of stress drops 
was shown to obey a power law 

Apart from comparing the statistics of stress drops from 
the model with that of the experimental time series, there 
is another more important issue, namely, does it generate 
power law statistics? If so, what mechanism is responsible 
for this? This is particularly important as the model is fully 
dynamical and noise free. It is clear that Fig|2ld is simi- 
lar to Fig. |2]b, as there no inherent scale in the magnitudes 
of the stress drops in both cases and thus it is likely to also 
show a power law statistics. Indeed, the distribution of stress 
drop magnitudes, D(A(|)), shown in Fig. |5]obtained from long 
runs for a large system size ( N — 1000) shows a power law 
Z)(A(|)) ^ A(|)^" over two orders of magnitude which increases 
with both length of stress series and system size. ( Note that 
the value of here nearly three times larger than the results in 
Ref. 1 32], Fig. 3b and thus, the power law is well converged 
with respect to the system size. ) Surprisingly, experimental 
points (•) corresponding to = 8.3 x lO^^i^' also fall on the 
same curve with an exponent value aw 1.1. ( We have scaled 
the experimental points by a constant amount along both the 
axis to show that these points also fall on the same line.) The 
distribution of the duration's of the stress drops D{At) ~ Af^P 
also shows a power law with an exponent value p w 1.3. The 
conditional average of A(|) denoted by < A(|) >c for a given 
value of Af behaves as < A(|) >c^ Af^/-^ with x w 0.65. The 
exponent values satisfy the scaling relation a = x(p — 1) + 1 
quite well. The exponent values remain unaltered in the region 
of strain rate 270 < e < 700 we have investigated thus is inde- 
pendent of the value of the drive parameter ( There are mod- 
els of coupled map lattices that produce power laws where 
the exponent value depends on the drive parameter.) We now 
investigate the underlying causes leading to this power law. 
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VI. DYNAMICS OF CROSSOVER 
A. Lyapunov Spectrum 

Our next aim is to characterize the dynamics of this 
crossover. A natural tool for characterizing the crossover is to 
study the distribution of Lyapunov exponents as a function of 
the applied strain rate in the entire interval where the PLC ef- 
fect is seen. Further, we also discuss the convergence proper- 
ties of the Lyapunov spectrum as the system size is increased. 
In particular, this will be useful in examining the density of 
null ( nearly vanishing) exponents and also to compare our 
results with that of the GOY model of turbulence. 

We have calculated the spectrum of Lyapunov exponents 
using the algorithm due to Benettin et al ll42il . The ex- 
ponent values reported here were obtained by averaging 
over 15000 time steps after stabilization with an accuracy 
of 10^*'. We have used several system sizes ranging from 

= 100, 150, 350, 500, 1000, and 3333 which covers approx- 
imately two orders of magnitude inM, i.e., from 301 to 10000. 
A rough idea of the changes in the dynamics of the system can 
be obtained by studying the dependence of the largest Lya- 
punov exponent (LLE) as a function of the strain rate. The 
LLE converges fast as a function of the system size. For in- 
stance, we find that the LLE for A^ — 500 looks much the same 
for a much smaller system size N=100 given in Fig. 3a in 
Ref. [32"]. The LLE becomes positive around e « 35 reaching 
a maximum at £ = 120, practically vanishing around 250. ( 
Periodic states are observed in the interval 10 < e < 35.) In 
the region 8 > 250, the dispersion in the value of the LLE is 
~ 5 X 10^"* which is the same order as the mean. Thus, the 
LLE can be taken to vanish beyond e = 250. 

The study of the Lyapunov spectrum reveals that in the 
chaotic regime of strain rates, only a small proportion of the 
exponents are positive, an equal small number are close to 
zero value and a large proportion of the exponents are neg- 
ative. The distribution of the Lyapunov exponents D(k), is 
shown for A^ = 1000 in the inset of Fig. |S]for the strain 
rate, e — 120. For this system size (with a total of exponents 
M=3001), the number of positive exponents is « 6.2% of the 
total number of exponents and the null exponents are also w 
9%. (For numerical purposes null exponents are taken to cor- 
respond to |A,| < 5.2 X lO^'*.) These ratios remain the same 
for the larger system sizes used. 

While D(k) reflects the distribution of Lyapunov expo- 
nents in various regions, for studying the convergence of the 
Lyapunov spectrum, plots of the density function j/Dky — 
/{XjOky /H) are better suited. Further, these quantities have 
been used traditionally in the studies of extended dynamical 
systems [271. A plot of J/Dky verses XjDky/H for e = 120 
for N = 500 and 1000 is shown in Fig. |6l It is clear that while 
the density function has not yet converged for negative values 
of XjDky/H, those for positive values are akeady converged. 
As we increase the strain rate beyond e = 180, concomitant 
with the decrease in the value of the LLE, the number of null 
exponents increases. For instance, at e = 220, for which the 
maximum Lyapunov exponent is small ^ 0.0058, the number 
of null exponents increases to 30% M (see inset of Fig. 0. 
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FIG. 6: A plot of j /Dky verses XjDky/H for e = 120. Inset shows 
a plot of D{X) as a function of X for N =1000. 



D{X) shows that the number of null exponents has increased. 
Concomitant with this trend, a plot of j/Dky verses XjDky /H 
for N = 500 and 1000 (Fig. Q shows that for e = 220 is weU 
converged for the entire range of values of the scaled Lya- 
punov exponent XjDky/H. This signals a faster convergence 
of the density function J/Dky = fO^jDKY /H) with the sys- 
tem size as we approach the scaling regime. Indeed, we find 
that plots for N = 500 and 1000 for strain rate e = 280 cannot 
be distinguished over the entire range of values of XjDky/H. 
Even though it would be adequate to use A^ = 1000, for further 
analysis when dealing with Lyapunov spectrum in the scaling 
regime, we use a much bigger system size of A^ = 3333, which 
for all practical purposes can be taken to be large A^ limit. A 
plot of J/Dky verses XjDky/H shown in Fig. |8]for A^ = 3333 ( 
and also for 1000) shows that for e — 280 is well converged for 
the entire range of values of XjDky /H. Note also that nearly 
40% the exponents are close to zero (see the inset). 

As we approach the power-law regime of stress drops (ex- 
tending from e = 250), as the largest Lyapunov exponent ap- 
proaches zero (~ 5.16 x 10^^ for e = 280), exponents below 
a certain value cross each other as a function of time. How- 
ever, the first few exponents remain distinct. Figure |9] shows 
the first two exponents that are well separated and another two 
which are close to each other in magnitude (for e ~ 280 and 
N = 3333). Below < 5 x 10"^ even though the expo- 
nents cross each other, the distribution of the exponents re- 
mains unchanged. The most significant feature of the spec- 
trum in the region is that there is a dense set of null exponents. 
The peaked nature of the distribution of the null exponents 
(|A.| < 5.2 X lO"'^) for e = 280 for A^ = 3333 is shown in Fig. 
M 

The peaked nature of for null exponents suggests the 

possibility of a power law distribution for their magnitudes. 
We have plotted the distribution of the null exponents (|A-, | < 
5.2 X 10^4) for e = 280, for a system size of A^ = 3333, M = 
10000 shown in Fig.[^ It is clear that both the positive and 
negative exponents show a power-law distribution £>(|A,|) ~ 
|A-|^^ with an exponent value y~ 0.51 and the scaling extends 
over an impressive three decades as shown in Fig. ^3 As null 
exponents correspond to marginal stable nature of the system. 
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FIG. 7: A plot of j/Dxr verses XjDKr/H for e = 220. Inset shows 
a plot of D{X) as a function of X for N =1000. 



FIG. 10: The peaked nature of the distribution of null exponents lying 
in the range [-5.2 x 10-'*,5.2x 10"'*] for £ = 280,^ = 3333. 
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FIG. 8: A plot of J/Dky verses XjDky /H for e = 280. Inset shows a 
plot of D(X) for N = 3333. A schematic plot of the Lyapunov density 
function (continuous line) for the GOY model (after ll28i1 V 




FIG. 1 1 : Log-log plot of the distribution of the marginal exponents 
for e = 280, N = 3333. Solid line is a guide to the eye. 



B. Comparison with Shell Model for IXirbulence 



their finite density, which itself obeys a power law, elucidates 
the underlying cause of power law distribution of stress drops 
at high strain rates. 
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FIG. 9: The first two Lyapunov exponents that do not cross each 
other as a function of time for N =3333 for e=280. Also shown are 
two more exponents that are close to each other. 



As mentioned in the introduction, both in the PLC ef- 
fect and in turbulence, power law statistics is seen at high 
drive rates in contrast to conventional SOC systems where the 
power law arises at low drives lEsil . In addition, the dense 
set of null Lyapunov exponents in the scaling regime is sim- 
ilar to the finite density of null exponents obtained by Ruelle 
1I43I1 for the discrete spectrum of the operator that linearizes 
the Navier-Stokes equations, if the fractal dimension of the 
energy dissipation set is Df < 2.5. This property is preserved 
by the GOY shell model |28]. Here we attempt a comparison 
of the Lyapunov spectrum obtained from our model with that 
of the GOY model. 

Shell models of turbulence fZT*] are designed to mimic the 
behavior of Navier-Stokes equations at high drives where the 
power law is seen. One standard model is the GOY model 
I27I l2^ . For this model, Ohkitani and Yamada fSs*] gave a 
good numerical evidence that the density function exists as 
the viscosity parameter r| tends to zero. In our case, the role 
of the viscosity parameter is taken by the applied strain rate. 
In Section IV A, we have shown that there is a rapid increase 
in the density of null exponents and consequently, there is a 
rapid convergence of j/Dky = fQ^jDKY /H) as function of 



N, starting from 8 = 220. This suggests that one should ex- 
pect convergence of the Hmiting j /D^y — fQ^jDKr/H) func- 
tion as we approach the power law strain rate regime of stress 
drops. Thus, we should expect that the limiting distribution 
itself converges as a function of e as we approach the scal- 
ing regime. Considering = 1000 approximates the limit- 
ing distribution ( see Fig. for justification), we have veri- 
fied that plots of J/Dky ~ fQ^jDKY/H) for three values of 
e = 250,260 and 280 for reasonably large N = 1000 con- 
verge. This result is similar to the convergence of the den- 
sity function in the GOY model as a function the viscosity 
parameter The density function j /Dky obtained from the 
model can be compared with that of the GOY model. Plot 
of i /Dky = .fO^jDKY /H) for a large system N=3333 ( which 
can be taken to represent the limiting density as a function 
of system size ) for e = 280 is shown in Fig. |8l along with a 
schematic plot for the GOY model shown by the continuous 
line. As can be seen, in both the cases, the distribution func- 
tion which is proportional to —df(k)/dX, shows a singularity 
near zero. The difference being that the singularity is more 
pronounced for our model. Ohkitani and Yamada also plot 
another quantity which represents the null exponents better, 
namely, the sum of Lyapunov exponents up to j normalized 
by H as function of j scaled by D^^y- The quantity J^l "^j/H 
is an increasing function of J/Dky for positive A-/ and goes 
to unity when ^^^j A,,- = H. In the region of null exponents, 
this quantity remains constant and then decreases with j when 
A-/s are negative. Thus, this quantity also reflects the density 
of null exponents. A schematic plot of ^ J '^j/H as a function 
of J/Dky ( continuous and dashed line) for the GOY model is 
shown in Fig. The increase in '^j/H for small j /Dky 
shows that there is a finite density of positive exponents in 
the Lyapunov spectrum for the GOY model. Further, these 
authors find that there is a convergence with respect to the 
viscosity parameter for Lyapunov spectrum corresponding to 
the interior of the attractor (ie., J/Dky < 1), while there is 
scatter for j /Dky > 1 (the dashed line represents this por- 
tion). We have plotted LjA,y/// as a function of j/Dky for 
N =3333 on the same plot for the sake of comparison. In our 
case, the increase to unit value is much slower ( compared to 
the GOY model) which clearly implies that there are very few 
positive exponents (of any significant magnitude) with most 
of them being vanishingly small. This feature is unlike the 
shell model where there is a finite density of positive expo- 
nents. In the GOY model, the largest exponent is proportional 
to r|^'/^ which is reflected in the steeper increase in Y,\ '^j/H 
for the GOY model. 

Thus, we have shown that the power-law regime seen in our 
model at high drives as in hydrodynamics is not a superficial 
feature. The Lyapunov spectrum of the model is quite similar 
to the GOY model which however has a finite density of pos- 
itive exponents. In addition, the distribution of the null Lya- 
punov exponents itself shows a power law in our case. Such a 
feature has not been verified for the GOY model. This feature 
is also quite distinct from the Lyapunov spectrum of models 
of SOC studied so far L44.i45i.i46i] . We comment on this later 
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FIG. 12: A plot of Y.{ /H as a function of J/Dky for e = 280 and 
A' = 3333. The corresponding schematic plot for the GOY model for 
1 1 Dky < 1 (continuous line) and J/Dky > 1 (dashed curve) (after 



C. Slow Manifold Analysis 

The analysis provided in the previous section shows that as 
the strain rate is increased most exponents get concentrated in 
a narrow interval around the zero value. As zero Lyapunov ex- 
ponents represent a marginal situation, we see that the region 
of strain rate beyond 240 (corresponding to the power law 
statistics of stress drops) can be identified with a marginally 
stable state. Thus, it would be interesting to realize a geomet- 
rical picture of dislocation configurations in the marginal state 
and examine how dislocations reach such state as a function 
of strain rate. We accomplish this through the slow manifold 
of the model O 

Recently, the geometry of the slow manifold of the origi- 
nal model has been analyzed in detail lITH IjtIi . The analy- 
sis shows that the relaxational nature of the PLC effect arises 
from the atypical bent nature of the manifold. Here we re- 
call some relevant results on the slow manifold of the original 
model {D — 0) and extend the ideas to the situation when the 
spatial degrees of freedom are switched on (D ^ 0). Slow 
manifold expresses the fast variable in terms of the slow vari- 
ables, conventionally done by setting the derivative of the fast 
variable to zero I13ll37ll 

p,n=g{p,„A) = -bopl+P,nd + pi,n=0. (9) 

where 5 = (|)"' — p,„, — a. The variable 5 has been shown 
to have all the features of an effective stress and thus plays 
an important physical role ^37^, particularly in studying the 
pinning-unpinning of dislocations. We note that 5 is a com- 
bination of two slow variables (|) and p„„ both of which take 
small positive values. Hence, 5 takes on small positive and 
negative values. Using Eq.|9] we get two solutions 

p,„ = [5+i8^+4boPi,ny^^]/2bo, (10) 

one for 5 < and another 5 > 0. For regions of 5 < 0, as bo 
is small ^ 10^'*, we get p,„/p,„, ~ —1/5 which takes on small 
values. This defines a part of the slow manifold, S2 where 
Pin is small. In this region, as the mobile density is small and 
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immobile density is large ( relative to p„,), this region can be 
identified with pinned configuration of dislocations and hence 
we shall refer to the region S2 as the 'pinned state of dislo- 
cations'. We note that larger negative values of 5 correspond 
to strongly pinned configurations, as they refer to smaller ra- 
tio of p,„/p„„. For positive values of 8, another connected 
piece is defined by large values of p„„ given by p„, w d/bo, 
which we refer to as the 'unpinned state of dislocations ' as 
Pirn is also small. These two pieces ^2 and are separated by 
5 = 0, which we refer to as the fold line ilStlSTiK see below). A 
plot of the slow manifold in the 5 — p,„ plane is shown in Fig. 
I13h . For the sake of illustration, we have plotted a monoperi- 
odic trajectory describing the changes in the densities during 
a loading-unloading cycle. The inset shows Pm{t) and 
For completeness, the corresponding plot of the slow mani- 
fold in the (p,„,p,„,,(|)) space is shown in Fig. 113b . along with 
the trajectory and the symbols. In this space, one can see that 
5 = (j)™ — p,„, — a = is a line that separates the pieces ^2 and 
Si of the slow manifold, and hence the name fold line. The 
cyclic changes in the variables is well captured by the nature 
of trajectory shown in Fig. 113b . The trajectory enters S2 at A 
and moves into S2, the value of 5 ( in Fig. 113b ) decreases from 
zero to a maximum negative value as the trajectory reaches B. 
Then 5 increases as the trajectory returns to A' before leaving 
52. The corresponding segment is ABA' in Fig. 113b . which is 
identified with the flat region of p„,(f) in the inset of Fig. 113b . 
As the trajectory crosses 5 = 0, dg/dpm becomes positive and 
it accelerates into the shaded region (Fig. I Oh ) rapidly till 
it reaches p„, = ?)/2bQ. Thereafter it settles down quickly on 

decreasing rapidly till it reenters 52 again at A. The burst 
in p„, (inset in Fig. 113b ) corresponds to the segment A'DA 
in Fig. 113b and b. The nature of trajectories for higher strain 
rate remains essentially the same, but is chaotic. The nature of 
the trajectories in the power law regime of strain rates which 
will be discussed later 

Having identified the regions of the slow manifold with the 
pinned and unpinned states of dislocations, we now consider 
the variation of stress when dislocations are pinned and are 
unpinned. Consider the stress changes as the state of the sys- 
tem goes though a burst of plastic activity. For D = 0, Eq. |6l 
reduces to 

(j) = <i[e-ep], (11) 

where = (j)'"p,„ defines the plastic strain rate. Since p,„ is 
small and nearly constant on ^2, stress increases monotoni- 
cally. However, during the burst in p,„ (A'DA in the inset), 
£,p{t) exceeds e leading to an yield drop. Since p,„ grows out- 
side S2, 5 = line separates the pinned state from the un- 
pinned state. Thus, 5 = physically corresponds to the value 
of the effective stress at which dislocations are unpinned. 

When the spatial degrees of freedom are included, there is 
no additional complication as the slow manifold is defined at 
each point. In this case, a convenient set of variables for vi- 
sualization of dislocations is (p,„(x),5(x),ji;). Here, our aim is 
to investigate the nature of typical spatial configurations in the 
chaotic and the power law regimes of stress drops and study 
the changes as we increase the strain rate. For simplicity, we 
shall use h — Q for which we have i^eff = (It is straightfor- 
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FIG. 13: (a)Bent slow manifold 5i and 52 (thick lines) with a simple 
trajectory for £ = 200 and m = 3. Inset: p,„ (dotted curve) and (|) 
(solid line).(b)The same trajectory in the ((]),p,m,Pm) space. 

ward to extend the arguments to the case when h ^ 0.) Then, 
the plastic strain rate tp(t) is given by 

ep(0=f"(f)7 /"'pm(x,fMx = r(f)p,„(f), (12) 
/ Jo 

where Pm(0 is the mean mobile density (= Y.j9m{jit)/N in 
the discretized form). With the inclusion of spatial degrees 
of freedom, the yield drop is controlled by the spatial average 
p,„(f) rather than by individual values of Pm{j)- Further, we 
note that the configuration of dislocations change during one 
loading-unloading cycle. However, one should expect that 
configurations will be representative for a given strain rate. 
Further, we know that the drastic changes occurs during an 
yield drop when p,„(f) grows rapidly. Thus, we focus our at- 
tention on the spatial configurations on the slow manifold at 
the onset and at the end of typical yield drops. 

First consider the configuration seen just before and after 
the yield drop when the strain rate is in the chaotic regime. 
In this regime, the stress drop magnitudes are large which im- 
plies that the change in mobile density is large. Figures[T4la. 
b for a typical value of e = 120. It is clear that both at the on- 
set and at the end of a typical large yield drop, the values 
which reflect the state of system ( pinned or unpinned state), is 
negative and correspondingly the mobile density p,„(j)'s are 
small, i.e., most dislocations are in a strongly pinned state. 



11 




-0.5 I -0.5 



FIG. 14: Dislocation configurations on the slow manifold at the inset 
and at the end of an yield drop: (a) and (b) for e = 120 (chaotic 
regime), and (c) and ( d) for e = 280 (scaling regime). 

( Recall that 5 signifies how close the spatial elements are 
close to unpinning threshold.) The arrows show the increase 
in Pm{j) at the end of the yield drop. We have checked that 
this is a general feature for all yield drops in the chaotic regime 
of strain rates. Now consider dislocation configuration in the 
scaling regime at high strain rates, say, e = 280, at the onset 
and at the end of an yield drop shown in Fig. Oc,d respec- 
tively. In contrast to the chaotic regime, in the scaling regime, 
most dislocations are clearly seen to he at the threshold of un- 
pinning with 5(j) w 0, both at the onset and end of the yield 
drop. This also implies that they remain close to this thresh- 
old all through the process of an stress drop. We have verified 
that the edge-of-unpinning picture is valid in the entire power 
law regime of stress drops for a range of values. Further, 
as a function strain rate, we find that the number of spatial el- 
ements reaching the threshold of unpinning 5 = during an 
yield drop increases as we approach the scaling regime. 

VII. TYPES OF BANDS 

The fact that the spatially extended Ananthakrishna's model 
is able to successfully reproduce the crossover dynamics from 
chaos to the power-law regime of stress drops ( and other 
generic features demonstrated earlier) might suggest that the 
characteristic features of the PLC bands may also emerge out 
of the model. Most models of dislocations bands use diffusive 
coupling although the physical mechanism of the term is dif- 
ferent in different situations 0] . An important feature of the 
spatial coupling in the model is that it accounts for spreading 
of dislocations into regions of low back stress once disloca- 
tions are unpinned ( the factor p,7„')- The term also determines 
the length scale over which dislocations spread into the neigh- 








FIG. 15: Spatially uncorrelated bands at e = 40. 



boring elements. Thus, while dislocation pinning and unpin- 
ning gives a heterogeneity in space (in principle), regions of 
low p,„, are favored for dislocation multiplication and spread- 
ing into neighboring regions. Further, this type of spatial term 
couples length scale and time scales in a dynamical way as 
p,m itself evolves in time and hence the associated time scale. 
Indeed, multiplication of dislocation depends on stress, {i.e., 
'^'"ff)' and hence this rate itself is changing dynamically lead- 
ing to changes in the time scale of internal relaxation as a 
function of e. We expect this to lead to changes in spatial 
correlation as strain rate is increased. 

Below we report both numerical and analytical studies on 
the spatiotemporal patterns emerging from the model as a 
function of the strain rate, e. We begin with the numerical 
results i47ll . 

For e < 10 and e > 2000, we get homogeneous steady state 
solutions for all the dislocation densities, p„„ p„„ and p^. In 
these ranges of strain rates, (|) takes the fixed point values 
asymptotically. In the region where interesting dynamics of 
chaotic and power law states are observed, the nature of the 
dislocation bands can be broadly classified into three differ- 
ent types occurring at low, intermediate and high strain rates 
described below. 

For strain rates, 30 < e < 70, we get uncorrelated static dis- 
location bands. The features of these bands are illustrated for 
a typical value, say for e = 40. A plot of p,„(j,f) is given in 
Fig- El Dislocation bands of finite width nucleate randomly 
in space and they remain static till another band is nucleated 
at another spatially uncorrelated site. The associated stress- 
time curves which are nearly regular have large characteristic 
stress drops. The distribution of these stress drops is found to 
be peaked as in experiments at low strain rates I20ll . 

At slightly higher values of strain rates, 70 < e < 180 we 
find that new bands nucleate ahead of the earlier ones, giving 
a visual impression of hopping bands. This can be clearly 
seen from Fig. 1161 where a plot of pm{j,t) is given for e = 
130. However, this hopping motion does not continue till the 
other boundary. They stop midway and another set of hopping 
bands reappear in the neighborhood. Often nucleation occurs 
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FIG. 16: Hopping type bands at e = 130 (arrow shows one such 
band). 



at more than one location. Stress-time plots in this regime 
have a form similar to Fig. ^ with the average amplitude of 
the stress drops being smaller than the localized non-hopping 
bands at low strain rates as seen in experiments. These stress 
drops also have a nearly symmetric peaked distribution as in 
the previous case but slightly skewed to the right similar to 
those observed in experiments |i,2Q.l . 

As the strain rate is increased further, the extent of propa- 
gation increases, concomitantly, the magnitudes of the stress 
drops decrease. We see continuously propagating bands even 
at e = 240 as can be seen from Fig. [n] One can see disloca- 
tion bands nucleating from one end of the sample (j = 0,t = 
25, 50 and 75) and propagating continuously to the other end. 
Often, we see a band nucleating at a point, branching out and 
propagating only partially towards both the ends. Unlike the 
present case which exhibits rather uniform values of p„,, we 
usually find irregularities as the band reaches the edges. The 
stress strain curves in this region of strain rates, exhibit scale 
free feature in the amplitude of the stress drops (Fig.|2]d) with 
a large number of small drops. As can be seen from Fig. |2]d, 
the mean stress level of these small amplitude stress drops in- 
creases until a large yield drop is seen. This large stress drop 
corresponds to bands having reached the end of the specimen. 

It is possible to calculate the velocity of the propagating 
bands in the high strain rate limit. We first note that our 
equations constitute a coupled set of integro-partial differen- 
tial equations, and hence cannot be dealt with in their present 
form. To reduce these equations to a form that is suitable for 
further analysis, we recall a few pertinent points about the 
changes in the structure of the slow manifold as a function 
of the applied strain rate. We note that the original model ex- 
hibits an incomplete approach to homoclinicity |37], i.e., the 
number of mixed mode oscillations of the type U, where s 
refers to the small period oscillations and L refers to a large 
relaxation oscillation, are limited. Typically, about 12 small 
period nearly harmonic oscillations are known to occur for a 
single large one at high values of the strain rate. The reason at- 








FIG. 17: Fully propagating bands at e = 240 
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FIG. 18: Slow manifold showing a trajectory for the space indepen- 
dent model near the reverse Hopf bifurcation point, at £ = 90, in = 2. 
• fixed point of Eqn 1-3 and 6. 



tributed to this is the finite rate of softening of the eigen value 
of the fixed point due to presence of the reverse Hopf bifurca- 
tion is reached IS 7,1 . In the presence of the spatial coupling we 
find that the softening is further enhanced as is clear from the 
fact that the upper Hopf bifurcation is pushed to much larger 
values of strain rate ( = 2000, see Fig. This enhanced 
softening rate implies that the number of small period oscil- 
lations is also increased in this domain of strain rates. Even 
so, the geometry of the slow manifold is not altered from that 
of the space independent model. In particular, the position 
of the unstable saddle focus remains located on the Si part 
of the manifold (Ref.ElIll). In addition, the feature of the 
fixed point approaching the fold line as a function of the strain 
rate is retained. Under these conditions, for high strain rates 
nearly sinusoidal oscillations are executed around the fixed 
point with the orbits touching ^2 only after executing several 
such turns. A plot of this is shown in Fig. ^] To understand 
the dynamics at high strain rates, we recall that our analysis 
in Ref. |37] shows that the orbit is re-injected along the sta- 
ble manifold close to the unstable saddle focus ( as shown in 
Fig. 11 ofRef.|37]). The orbit then spirals out along the unsta- 
ble manifold of the fixed point. Once the orbit is sufficiently 
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away from the fixed point when the influence of the fixed point 
is lost, it is re-injected close to the fixed point via Si. The dy- 
namics then repeats. Note that at high applied strain rate, the 
system is close to reverse Hopf bifurcation point and hence the 
fixed point is close to the applied strain rate value. Thus as the 
orbit executes one turn, there is one small yield drop. How- 
ever, the orbit executes several turns around the fixed point, 
each turn leading larger loop sizes, i.e., larger values of p„, 
and consequently to successively larger stress levels than the 
earlier one before briefly visiting ^2. 

Under these conditions the dynamics is entirely controlled 
by the spiralling motion around the fixed point. Thus, the en- 
tire dynamics is essentially described by the fast variable; the 
other two variables p„„ and (|) can be taken to be parameters. 
Such a situation is described by the transient dynamics dic- 
tated entirely byequation of the fast variable (the so called 
layer problem IBSll ') and thus, we are justified in using only 
the evolution equation of the fast variable in terms of the slow 
manifold parameter 5 = (|)"' — Pim — a. Since the trajectory 
rarely visits the ^2 part of the slow manifold, we restrict the 
calculations to 5 > 0. The physical picture of a propagating 
solution is that as the orbit at a site makes one turn around 
the fixed point, i.e., 5 small but positive, around the value of 
the applied strain rate, the front advances by a certain distance 
along the specimen like the motion of a screw. 

The rate equation for the mobile dislocation density p„, in 
terms of 5 is 



9pm 

dt 



oVm 



-boP 



-?>Pm + Pim+D' 



, d^Pm 

3x2 



(13) 



where D' = D(|)"Yp„„. Since, the slow variables, p,„, and (|) are 
treated as parameters, this has the form of Fisher- Kolmogorov 
equation for propagating fronts, which has been well studied. 
This equation can be reduced to the standard form 



dZ 



--Z{l-Z)+D 
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(This is done by first transforming p,„ ~X — p,„i/5, dropping 
the term IboPi^/b compared to 5 in the linear term in X, and 
then using Z = Xb/bo and t' = f5.) It is clear that Z = is 
unstable and Z = 1 is stable. Using the form for propagating 
front Z = Z„e™'-*-'^', the marginal velocity is calculated us- 
ing V* — Re(i>{k*) / Rek* = d(i>/ dk\i^=i^* and 3t/co/£/^|j:=t:* =0, 
gives the velocity of the bands v* = 2 f4^ l49il . In terms of the 
variables in Eqn.[OJ the marginal velocity is 
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In order to relate this to the applied strain rate, we note that for 
a fixed value of the strain rate (where propagating bands are 
seen), the average level of stress drop is essentially constant. 
Thus, from Eqn. |6l we see that in this regime of high strain 
rates, the applied strain rate 8 is essentially balanced by the 
plastic strain rate (l//) /o(|)"'p,„(x,r) = tp. Then, using (|)™ = 
e/pm, and using 5 = (|)"' — a — p,„„ we get 



I Dt , t 
{- a -pi, 

PmPim Pm 



It is important to note that at high applied strain rate p^ p*, 
, the fixed point value. Thus, for all practical purposes, we 
can assume p,„ as a constant. From the above equation, we 
see that the velocity of the propagating bands is proportional 
to the applied strain rate. This result is similar to the result 
obtained recently by Hahner et al. llsoll . Further, v ^ p^^ 
which also appears to be consistent with an old experimental 
result. (See Fig. 7 of Ref. |51] which appears to fit vp„, = 
constant.) This result needs further experimental support. 

As the form of our equation has been reduced to the stan- 
dard form, all other results carry through, including nonlin- 
ear analysis. We have numerically calculated the velocity of 
the continuously propagating bands at high strain rates from 
the model which confirms the linear dependence of the band 
velocity on applied strain rate. In the region of strain rates 
e = 220 to 280 (corresponding to unsealed strain rate values 
10-5- 1.5 X lO^V^), we find that the unsealed values of the 
band velocity increases from 100 to 130 pm/s. These val- 
ues are consistent with the experimental values reported by 
Hahner et. al. jsoll. 
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We note here that the types of the bands seen in our model 
are correlated with the two distinct dynamical regimes inves- 
tigated. The hopping type bands belong to the chaotic regime, 
a result consistent with the recent studies on Cu-Al polycrys- 
tals 1 20]. On the other hand, the propagating bands are seen 
in the power law regime of stress drops 1 32], again consistent 
with these studies ] 19, 20]. Curiously the uncorrected bands 
predicted by the model also belong to the chaotic regime. We 
shall now explain these results based on the dynamics of the 
model. We first note that each spatial element is described 
the three dislocation densities ( Eq. (1-3)). Consider one of 
these elements being close to unpinning threshold, ie., 5 = 0. 
It has been shown earlier that p„„ is out of phase with p,„ 
lITsl Wh . This feature is retained with the spatial coupling 
as well. When the orbit is about to leave ^2, ie., when Pm{j) 
is at the verge of a sharp increase, p„„ is largest. However, 
the extent of the spatial coupling is determined by pj^- But 
the magnitude of p,„, itself decreases with the applied strain 
rate, being large at low strain rates 113 . 37.1 . Thus, the spatial 
width of this is small at low e and large at high e. Next we 
note that the growth and decay of Pm{i) with j occurs over 
a short time scale which is typically of the order of the cor- 
relation time, ic, of (|)(f). Beyond this time, the memory of 
its initial state is lost. Consider an initial state when a band 
is formed at some location. Before the memory of this initial 
state decays, if a new band is not created, we get an uncorre- 
cted band. On the other hand, if a new band is created before 
the memory of the initial state decays, there are two possi- 
bilities. If another band is created just before the correlation 
decays substantially by that time, we get a hopping type band. 
If however, even before the burst of p,„(y) decreases beyond 
its peak value, new sources of creation of p,„ occur, then we 
end up seeing a propagating band. An analysis of the corre- 
lation time shows that it increases with the applied strain rate. 
Concomitantly, p„„ decreases with e which implies that the 
spatial correlation increases. (Indeed, the value p,„, is quite 
small for large £ as we reach the power law regime of stress 
drops.) Under these conditions, only partial plastic relaxation 
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is possible in this regime. This discussion clarifies the dy- 
namic interplay of time scales and length scales. Moreover, as 
the spatial coupling term allows the spreading of dislocations 
only into regions of low p,„, or low back stress, the propen- 
sity for continuous propagation of the band is enhanced when 
Pi,„ is small. In addition, we find that higher values of p„„ at 
the wake of the band which favors propagation into regions 
of smaller immobile density thus determining the direction of 
propagation also. 



VIII. SUMMARY AND CONCLUSIONS 

We first summarize and make appropriate comments wher- 
ever necessary. Detailed numerical and analytical studies on 
the extended Ananthakrishna's model shows that it repro- 
duces all the important features of the PLC effect including 
the crossover in the dynamics from a chaotic to a power law 
regime observed in experiments. It also provides insight into 
the dynamical causes leading to this crossover. A systematic 
study of the system size effects of the Lyapunov spectrum car- 
ried out elucidates the underlying mechanism controlling the 
crossover. The study demonstrates that the limiting Lyapunov 
distribution evolves from a set of positive and negative expo- 
nents with a few null exponents to a dense set of null expo- 
nents as we approach the scaling regime of stress drops. This 
study is complemented through an analysis of the slow man- 
ifold. This method is particularly useful in giving a geomet- 
rical picture of the spatial configurations, both in the chaotic 
and power law regime of stress drops. The study shows that 
the configuration of dislocations is largely in the pinned state 
in low and medium strain rates (chaotic domain) are pushed 
to the threshold of unpinning as we increase strain rate ( 
power law stress drop regime). The study also establishes 
that the present model has considerable similarities with the 
GOY model of turbulence [28 1. The model also reproduces 
the major spatial features of the PLC effect. The randomly 
nucleated band, the hopping and propagating types are found 
as the strain rate is increased. It also predicts a linear depen- 
dence of the velocity of the band and inverse dependence on 
the mobile density at high strain rates. 

Several observations may be in order on the dynamics of 
the crossover. We first note that the crossover itself is smooth 
as the changes in the Lyapunov spectrum are gradual, though 
it occurs in a narrow interval of strain rates ( from 220-250). 
Second, the power law here is of purely dynamical origin (in 
the sense elaborated below below). We have shown that this 
is a direct result of the existence of a reverse Hopf bifurca- 
tion at high strain rates. In this regime due to softening of the 
eigen values ( as a function of the applied strain rate), the or- 
bits are mostly restricted to the region around the saddle node 
fixed point located on the Si part of the manifold. This offers 
a dynamical reason for the smallness of the yield drops in this 
region fTslls^ . Note also that there is a dynamic feed back 
between the stress determined by Eq. |6l and the production 
of dislocations in Eq. ^ which provides an explanation for 
the slowing down of the plastic relaxation. ( The partial plas- 
tic relaxation has been cited as the reason for the power law 



.) This sets up a competition between the time scale 
of internal relaxation and the time scale determined by the ap- 
plied strain rate ( essentially Deborah number). We note that 
while the time scale for internal relaxation is increasing, that 
due to the applied strain rate is decreasing. Third, our analysis 
shows that the power law regime of stress drops occurring at 
high strain rates belongs to a different universality class com- 
pared to SOC systems, as it is characterized by a dense set of 
null exponents. This must be contrasted with the lack of any 
characteristic feature of the nature of the Lyapunov spectrum 
in the few models of SOC studied so far |44,l45l l4al . For in- 
stance, no zero and positive exponents, zero exponent in the 
large limit etc., have been reported |44j, j45|, ,46]. ( Often, 
the nature of largest Lyapunov exponent is inferred based on 
the similarity of other dynamical invariants f^T].) The dense 
set of null exponents in our model is actually similar to that 
obtained in shell models of turbulence where the power law 
is seen at high drive values |28]. However, there are signifi- 
cant differences. First, we note that the shell model |28| can- 
not explain the crossover as it is only designed to explain the 
power law regime. Second, the maximum Lyapunov exponent 
is large for small viscosity parameterrj, ie., A-i °^ rj^'/^ in shell 
models |28] in contrast to near zero value in our model. It is 
also interesting to note that in our model propagating solu- 
tions arise in the power law regime of stress which comes as a 
surprise. As far as we are aware, this is the first situation, both 
from a experimental and theoretical angle, where propagating 
solutions are seen in a marginally stable situation. 

Regarding the spatial features seen in the model, we stress 
that these features emerge purely due to dynamical reasons 
without any recourse to using the negative strain rate sen- 
sitivity feature as an input, as is the case in most models 
L50, 52. 53. ,54,1 . Even the recently introduced poly-crystalline 
plasticity model which reproduces the crossover behavior also 
uses the negative SRS as an input fsS*]. The dynamical ap- 
proach followed here clearly exposes how the slowing down 
of the plastic relaxation occurs due to a feed back mecha- 
nism of dislocation multiplication and applied strain rate as 
we reach the power law regime of stress drops. While the 
three different types of bands have features of the uncorre- 
cted type C, hopping type B and the propagating type A bands 
found in poly=crystalline materials, there is no element of 
poly-crystallinity in the model in its present form. In poly- 
crystals, other types of coupling terms do arise which have 
also been modelled by diffusive type terms (l|. One way of 
including the presence of grain boundaries within the natu- 
ral setting of the model is to recognize that crossslip will be 
hindered near the grain boundaries which also leads to a term 
similar to the present diffusive term. Such a term can account 
for the back stress arising from the incompatibility of grains. 
As the form of these terms are similar, the basic results are 
unlikely to change although one should expect a competition 
between the terms operating within a crystal and that at the 
grain boundaries. 

From a purely dynamical point of view, this model should 
be of interest to the area of dynamical systems, as it appears to 
be the first fully dynamical model which exhibits a crossover 
from a chaotic to a power law regime, in the sense that our 
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model is continuous space time model without any recourse 
to artificial thresholds as is done in coupled map lattices is^l . 
We note also that while the slow manifold subspace gives a 
method of visualizing the dislocations configurations, partic- 
ularly in the scaling regime, the complementary subspace of 
the fast variable has helped us to obtain the band velocity in 
the same regime of strain rates. From the point of view of 
plastic instabilities, the present dynamical approach should be 
a promising direction for explaining many other patterns men- 
tioned in the introduction frl. 

Finally, as stated in the introduction, the PLC effect bears 
considerable similarity with many stick-slip systems and 
hence presents a way of understanding some of these systems. 
Here we make a few comments on the similarity of the present 
model for a possible adoption to the observed voltage fluctua- 
tions in charge density wave (CDW) compounds 1 11, 591. Un- 
der the action of applied electric field, anomalously large volt- 
age fluctuations are reported when the electric field is above 
the threshold value. This ohmic to non-ohmic transition in 
Kq^MoOj, and Rb(j T,Mo03, for instance |57], has not been ad- 
equately explained although the similarity with the PLC effect 



has been noted fT]\. Lee and Rice f5§\ have suggested that 
phase dislocations of the CDW carry current at fields too low 
for the CDW to move as a whole. Indeed both fall in the cate- 
gory of pinning-depinning phenomenon. In the case of CDW, 
pinned at impurity/defect sites is unpinned due to the applied 
electric field. The threshold value of the electric field can be 
viewed as the onset of plastic flow of the charged phase dis- 
locations 1 11]. These authors identify stress with voltage and 
strain with current and suggest that the total current is the sum 
of the ohmic part and that arising from charge density waves. 
This is corresponds to the elastic and plastic displacements 
competing to give rise to the PLC instability. We believe that 
this parallel can be taken further along the lines of our PLC 
model where one can identify the phase dislocations with mo- 
bile dislocations, the neutral defects of the CDW with dislo- 
cation dipoles, ie., the immobile, and the phase dislocations 
pinned at defects with the Cottrell type. Work along these 
lines is progress. 
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